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ABSTRACT 

We discuss the properties of the X-ray sources with faint optical counterparts in the 
very young open cluster NGC 6231. From their positions in the H-R diagram, we 
find that the bulk of these objects probably consists of low-mass pre-main sequence 
stars with masses in the range 0.3 to 3.0 M . The age distribution of these objects 
indicates that low-mass star formation in NGC 6231 started more than lOMyr ago and 
culminated in a starburst-like event about 1 to 4Myr ago when the bulk of the low- 
mass PMS stars as well as the massive cluster members formed. We find no evidence 
for a spatial age gradient that could point towards a sequential star formation process. 
Only a few X-ray sources have counterparts with a reddening exceeding the average 
value of the cluster or with infrared colours indicating the presence of a moderate near- 
IR excess. The X-ray spectra of the brightest PMS sources are best fitted by rather 
hard thermal plasma models and a significant fraction of these sources display flares 
in their light curve. The X-ray brightest flaring sources have decay times between 2 
and 16 ks. The X-ray selected PMS stars in NGC 6231 have logLx/^boi values that 
increase strongly with decreasing bolometric luminosity and can reach a saturation 
level (log Lx/L\y Q i ~ —2.4) for non-flaring sources and even more extreme values during 
flares. 

Key words: Stars: pre-main sequence - X-rays: individuals: NGC 6231 - X-rays: 
stars - Open clusters and associations: individual: NGC 6231 



1 INTRODUCTION 

The formation of massive stars is currently one of the key 
questions in stellar astrophysics. Various scenarios have been 
proposed; accretion from a circ umste llar envelope (e.g. Shu 
et al. Il987l , Behrend & Maeder 1200 ll ) , accre tion from a tur- 
bulent molecular core (McKee & Tan l2003l) or a mixture of 
competitive accretion and collisions of lower mass protostars 
in the core of de nse st ellar clusters (e.g. Bonnell et al. 1 19981 . 
Bonnell & Bate 120021 ) being the most popular ones. Since 
the bulk of the massive stars are found in open clusters, 
some clues on this question can probably be obtained from 
a study of the star formation processes within young open 
clusters. Beside their massive star population, such clusters 
usually harbour a wealth of lower mass objects with the 
least massive ones that have not yet reached the zero age 
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main sequence. An important issue is the relationship be- 
tween the high-mass cluster members and these low-mass 
objects. While it has sometimes been suggested that the 
massive stars trigger the formation of the lower mass ob- 
jects, in several cases, the lower mass pre-main sequence 
stars are found to be significantly older than the massive 
clust er membe rs (e.g . Damiani et al. |2004| . Prisinzano et al. 
120051) . Herbig j 19621 ) accordingly suggested that the forma- 
tion of early-type stars could stop the formation process of 
lower mass stars in a cluster by disp ersing the ambient gas. 
Other authors (e.g. Sung et al. 1 19971 ) however suggested that 
the observed age discrepancy might result from the fact that 
the PMS evolutionary tracks actually overestimate the real 
age of the protostars. Obviously, investigating the proper- 
ties of lower mass pre-main sequence (PMS) stars in very 
young open clusters can help us to understand the feedback 
of massive stars on their environment and can thereby shed 
light on the way the most massive objects have formed. 

Classical TTauri PMS stars display emission in the 
Balmer lines and can thus be identified through photometric 
observations with an Ha narrow band filter. Most classical 
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T Tauri stars (cTTs) also display near-infrared excesses that 
are attributed to heat ed du st in a disk-like accretion struc- 
ture (e.g. Meyer et al. Il997h . However, another category of 
PMS objects, the so-called weak-line T Tauri stars (wTTs), 
display only weak optical line emission and do not show 
significant near-IR excesses. An interesting property that 
can help us to identify these PMS stars is their relatively 
strong X-ray emission with l og Lx /Lboi reaching values as 
large as —3 (e.g. Neuhauser 1 19971 . Feigelson & Montmerle 
1 19991 ). X-ray observations can therefore enable the selection 
of PMS stars, such as wTTs, that would neither be detected 
throu gh Ha nor through near-IR photometry (e.g. Damiani 
et al. l2004h . The high sensitivity of XMM-Newton and the 
exceptional spatial resolution of Chandra have already been 
used to study some rich young open clusters in the X-ray 
domain. For instance, NGC 6530 at the core of the Lagoon 
Nebula was found to harbour both classical and weak-line 
T Tauri stars with masses between 0.5 and 2.0 M®, but only 
a few of the X-ray s electe d PMS candida tes we re found to be 
cTTs (Rauw et al. |2002| . Damiani et al. l2004h FI NGC 6383, 
on the other hand, does not contain any known cTTs, but 
our XMM-Newton observations revealed a number of weak 
X-ray sources in this clus ter th at are most likely associated 
with wTTs (Rauw et al. 120031 ). A preliminary analysis of 
the optical properties of these sources suggested that they 
are in fact older than the massive binary HD 159176 (07 V 
+ 07 V) in the centre of the cluster. 

NGC 6231 - at the core of the Sco OBI association 
- is another very young open cluster, rich in massive O- 
type stars. Using UBV(I)cHa photometry, Sung et al. 
(j 19981 . hereafter SBL98) inferred a mean cluster redden- 
ing of E(B - V) = 0.466 ± 0.054 and a distance modu- 
lus of 11.0 ± 0.07. These authors noted that the redden- 
ing law towards NGC 6231 could be somewhat peculiar with 
Rv = 2.45 E(V - Ic)/E(B - V) = 3.3 ±0.1. SBL98 found 
only 12 PMS objects (and 7 PMS candidates) brighter than 
V = 17 displaying Ha emission. However, as pointed out 
above, there could be a number of PMS stars without Ha 
emission; but with optical photometry only, these would be 
very difficult to distinguish from field stars lying to the right 
of the cluster main-sequence. In this paper, we use the cat- 
alogue of XMM-Newton sources in NGC 6231 presented by 
Sana et al. (2006a], hereafter Paper I) to investigate the prop- 
erties of X-ray sources with faint optical counterparts in 
this cluster. The X-ray properties of the early-type stars in 
NGC 6 231 were discussed in a separate paper (Sana et al. 
l2006bl . Paper II). 



2 OPTICAL COUNTERPARTS 

Adopting a 3" correlation radius, 446 of the X-ray sources 
in the NGC 6231 field h ave an optical counterpart in the 
catalogue of Sung et al. (1 19981 ) extended down to V = 21 

We note here that the luminosities derived by Rauw et al. 
(20Q2) for the PMS objects in NGC 6530 are too low, leading 
to an overestimate of their age. The reason for t his er ror is that 
the V — Ic colour index quoted by Sung et al. ll2000h are actu- 
ally given in the Cousins photometric system while the Vizier 
database erroneously indicates that these are expressed in the 
Johnson system. 



(Sung H., unpublished). We note that our combined EPIC 
field of view extends over a wider area than the actual field 
investigated by SBL98 and some sources outside the SBL98 
area actually have counterparts in other optical/near-IR cat- 
alogues (see Paper I). Nevertheless, throughout this paper, 
we will restrict ourselves to the sources with counterparts 
in the extended SBL98 catalogue. The reasons are, on the 
one hand, that the USNO and GSC photometries are less 
accurate than the results of SBL98. On the other hand, the 
new photometry acquired by one of us (H. Sung) over a 
larger field than initially investigated by SBL98 (referred to 
as SSB06 in Paper I) does not include observations with the 
Ha filter. 

The V vs. B — V and V vs. V — Ic colour-magnitude 
diagrams of the optical counterparts are shown in Fig.[T] 
The bulk of the objects are found to be rather faint stars 
mostly located to the right of the main-sequence. At V — 18, 
the photometric error on most colour indices reaches about 
0.1 mag. Errors on the B — V index tend however to be 
larger (Fig. [1] upper panel) because of the lower quality 
of the B photometry in SBL98. This fortunately does not 
affect our work as we focus on the V, V — Ic and R — Ha 
measurements. 

Among the 446 counterparts in the extended version of 
the SBL98 catalogue, 392 sources have been observed in the 
V , Ic and Ha bands. This allowed the authors to compute 
a pseudo continuum (labelled R) near Ha, as the mean of 
the V and 7c magnitudes. To identify Ha emitting stars 
among these X-ray selected objects, one can then compare 
the R — Ha index to its value for main sequence stars. For 
the latter, we use the relation between the R — Ha and 

V — Ic in dices for main sequence stars proposed by Sung et 
al. (|l997h and we account for the effect of cluster redden- 
ing on the V — Ic colour index (see below). Following Sung 
et al. (|l997T ). we consider that a star displays Ha emission 
if A(R - Ha) = (R - Ha) - (R - Ha) M s 5* 0.21 mag, 
whereas a star will be considered an Ha emission candidate 
if 0.12 < A(R - Ha) < 0.21. 19 objects of our sample have 

V — I < 0.50 and one has V — I > 3.0 . Hence, they do 
not fall in the range covered by the (R — Hu)ms relation 
of Sung et al. (ll997T l. From a visual interpolation (see Fig. 
2, upper panel), one can still conclude that these are defi- 
nitely not Ha emitters. In this way, we find that about 18% 
(70 out of 392) of the X-ray selected objects are either con- 
firmed or potential Ha emission objects. This large number 
of Ha emission candidates compared to the original results 
of SBL98 is mainly due to the fact that the photometry used 
here extends down to much fainter magnitudes (Fig. [2] lower 
panel). 

A total of 384 EPIC sources in the field of view of 
NGC 6231 have also a counterpart in the Two Micron All 
Sky Survey point source catalogue (2MASS, Cutri et al. 
2003). Among these, 333 sources have also a counterpart 
in the SBL98 catalogue. We used the March 2003 update 
of t he col our transformations, initially derived by Carpen- 
ter (|200lT) and available on the 2MASS websitcQ, to convert 
the J — H and H — K s colour indices to the homogenized 
JHK photometric system introduced by Bessell & Brett 
(1988). Assuming that the 2MASS and SBL98 magnitudes 
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Figure 1. Colour-magnitude diagrams of the X-ray sources 
with an optical counterpart in the SBL98 catalogue. Filled dots, 
filled diamonds and open triangles indicate respectively Ha emit- 
ting stars, Ha candidates and stars with no evidence for emis- 
sion (see Fig.[2]|- The reddening vector wit h Ry = 3.3 and 
E(V - I C )/E(B - V) = 1.365 (Sung et al. Il99gft is indicated 
and the solid line shows the ZAMS relation taken from Sung & 
Bessell shifted by a distance modulus DM = 11.07 and 

reddened with E(B - V) = 0.466 (Sung et al. ll998h . 

are not affected by photometric variability, we tried to estab- 
lish the effective temperature and bolometric correction us- 
ing a dereddened V—K colour inde>0. Comparing the results 
with those obtained from the V — Ic index (see below), we 
find that, for the large majority of the stars, the former tech- 
nique yields systematically lower temperatures (sometimes 

3 For this purpose, we use the reddening law of Rieke & Lebofsky 
lll985ll . 
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Figure 2. Top panel: R — Ha as a function of V — Ic for 
the SBL98 optical counterparts. The solid line yields the R — Ha 
ver sus V — Ic relation for main sequence stars taken from Sung et 
al. lll997f) and reddened with the average reddening of NGC 6231. 
The short- and long-dashed lines yield respectively the thresholds 
for Ha emission candidates and Ha emitters (see text). Bottom 
panel: V-magnitude as a function of R — Ha for the SBL98 
optical counterparts of our EPIC sources. The symbols have the 
same meaning as in Fig.[T] 

by more than 1000 K) than the latter one. For objects sur- 
rounded by circumstellar material, such a discrepancy could 
reflect the presence of a near-IR excess (that would mimic a 
lower temperature in the V — K colour index) . However, for 
the bulk of the objects, the most likely explanation for this 
effect is probably the rather large errors in the 2MASS mag- 
nitudes. Note also that the reddening correction over such 
a wide wavel ength range, using the extinction law of Rieke 
& Lebofsky (|l985l ). might simply not be adequate for stars 
in NGC 6231 (perhaps as a result of the peculiarities of the 
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Figure 3. JHK colour-colour diagram of the 2MASS counter- 
parts, with good quality near-IR photometry, of the X-ray sources 
in the EPIC field of view around NGC 6231. The heavy solid line 
yields the intrin sic colours of main sequence stars following Bessell 
& Brett l l 1988ft . whereas the reddening vector is illustrated for 
Ay = 1.538 and adopting the reddening law of Rieke & Lebofsky 
l|l985h . The crosses stand for 2MASS sources having no counter- 
part in the SBL98 catalogue (mainly because they fall outside the 
field of view investigated by Sung et al. I1998T ). The other sym- 
bols have the same meaning as in Fig.[TJ The dotted straight line 
yields the locus of dereddened colours of classical TTauri stars 
according to Meyer et al. 



extinction law discussed by Sung et al. 1 19981 ). In addition, 
the 2MASS colour indices and magnitudes of some objects 
are either subject to large uncertainties or are only upper 
limits due to non-detections. In total, we find that only 295 
out of the 384 2MASS counterparts have quality flags A, B, 
C or D for the measurements of all three individual near-IR 
magnitudes. Restricting the JHK colour-colour diagram to 
these objects, we find that their location is consistent with 
slightly reddened main-sequence or giant stars (see Fig. [3}. 
Only a couple of sources show evidence for a moderate IR 
excess and only about ten objects display strongly reddened 
IR colours. 

Although a small fraction (~ 9%, see Paper I) of the 
X-ray selected objects might in fact be field stars (either 
foreground or background) unrelated to NGC 6231, we have 
assumed that all stars are located at the distance of the 
cluster (DM = 11.07 ±0.04, Paper I) and a re all subject to 
the same reddening (Av = 1.538, Sung et al. 1998). We have 
then built the Hertzsprung-Russell diagram (HRD) of these 
EPIC sources using effective temperatures and bolometric 
corrections interpolated from the dereddened V — Ic colour 
indices for main sequence stars of s pectr al types B0 to M6 
tabulated by Kenyon & Hartmann (1995). Eight sources (a 
subset of those not covered by the [1997] Sung et al. relation) 
do not fall in this range and are therefore not plotted in the 
HRD. 

Figure[4] compares the location of the X-ray selected ob- 



Figure 4. Hertzsprung-Russell diagram of the EPIC sources with 
optical counterparts in the SBL98 catalogue. The symbols have 
the same meaning as in Fig. CO Evolutionary tracks from Siess et 
al. ll2000h for masses of 0.2, 0.3, 0.4, 0.5, 0.7, 1.0, 1.5, 2.0, 2.5, 3.0, 
4.0, 5.0, 6.0 and 7.0 Mq are overplotted. The thick solid line shows 
the ZAMS. The thick dashed line indicate the birth isochronc 
(i.e. the Myr isochrone) from Siess et al. while the other dashed 
lines correspond to isochrones for ages of 0.5, 1.5, 4.0, 10.0 and 
20.0 Myr. The stars located to the left of the vertical dashed line 
were not considered in the age distribution of the cluster. 



jects t o the pre-main sequence evolutionary tracks of Siess 
et al. (|2000|) for Z = 0.02 and without overshooting (note 
that these evolutionary models include neither rotation nor 
accretion). Errors on the effective temperatures and bolo- 
metric magnitudes were obtained by error propagation, in- 
cluding contributions from the photometric measurements, 
the cluster distance, and the average E(B — V) and E(V — I) 
color excesses obtained from SBL98. Note that the color de- 
generacy of the hotter stars (i.e. the fact that a small range 
in color index corresponds to a large range of T e ff) implies 
that the determination of their effective temperature be- 
comes very sensitive to any uncertainty on the (V — J)o 
index. Even when the observed V — I data are excellent, the 
dereddening introduces a fix contribution to the error on 
(V — I)o, yielding uncertainties on the effective temperature 
of several thousands Kelvin (see Fig. [4|. In the following, 
the twelve stars with T e ff > 12500 K will be ignored while 
computing the age distribution of the stars in NGC 6231 as 
they do not provide any constraint on the cluster age. 

We find that most of the X-ray selected, optically faint 
objects (V ^ 14) fall between the PMS evolutionary tracks 
for stars of masses in the range 0.3 to 3.0 Mq. A compari- 
son with the isochrones further indicates that the majority 
of these objects should have ages ranging from 1 to 10 Myr 
(see Fig.[5]). We estimated the age of the 372 objects by lin- 
early interpolating the isochrones (computed between and 
20 Myr with a step of lMyr) along the M bo i axis. The ob- 
tained distribution of ages peaks around ~ 1 - 4 Myr and 
decreases towards lower ages (see Fig. [5}. Average and quar- 
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Table 1. Average, dispersion and quartiles (in Myr) of the age 
distributions of the PMS populations in NGC 6231. N indicates 
the number of objects considered in each (sub-)population. 



Population 


N 


mean 


1-cr 




Quartiles 








dispers. 


1st 


2nd 3rd 


wTTs 


302 


4.6 


3.7 


1.9 


3.4 4.9 


cTTs 


70 


4.8 


1.5 


1.9 


2.9 3.3 


All 


372 


4.6 


3.9 


1.9 


3.5 5.6 



tiles are given in Table [T] The mean and 1-a dispersion are 
strongly affected by the older age tail and does not represent 
neither the centroid nor the width of the peak seen in the 
histograms. Within our uncertainties, no significant differ- 
ence is found between the different distributions considered 
in Tabled] 

We note that Sung et al. (119981 ) compared the H-R di- 
agram of NGC 6231 to PMS tracks of Bernasconi & Maeder 
(1996) finding an age spread of about 11 Myr for the low 
mass objects. While this is similar to the age spread found 
here, we caution however that a direct comparison of age 
determinations of PMS stars made with different evolution- 
ary models and different calibrations might be biased (see 
e.g. the discussion in Siess et al. I2000T ) . As for NGC 6530, 
we do not find a large age difference between those stars 
displaying Ha emission (probable cTTs) and those with- 
out (probable wTTs). What we do see, is that Ha emission 
seems restricted to X-ray selected stars with masses below 
about 2.5 Mq. This is again reminiscent of the situation in 
NGC 6530. 

Zinnecker (120021 ) noted that the binary frequency 
among pre-main sequence stars is at least as high as among 
main sequence stars. Binarity affects the distribution of 
a coeval population of PMS stars in the H-R diagram by 
introducing a ban d shifted upwards from the true isochrone 
(Siess et al Il997l ) leading to an underestimate of the actual 
age. Palla (2002) evaluated the average age discrepancy for 
a realistic distribution of binary mass ratios, finding that 
age estimates that do not account for binarity can be off 
by a factor ~ 1.5. Therefore, binarity should not have too 
large an impact on our age estimates above. In any case, 
binarity cannot explain the tail of the age distributions 
towards older ages (see Fig. [5}. 



3 X-RAY PROPERTIES OF THE FAINT 
SOURCES 

3.1 X-ray spectra 

Spectra were extracted for each source, for each instrument 
and each observation using the procedure described in Pa- 
per II. Briefly, we adopted circular extraction regions with 
a radius corresponding to half the distance to the near- 
est neighbouring X-ray source. For each source, the back- 
ground spectrum was obtained over source-free regions cho- 
sen according to the source location. We adopted the redis- 
tribution matrix files provided by the XMM-Newton SOC, 
whereas the ancillary response files were built with the ap- 
propriate SAS tasks. The spectra were binned to have at 
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Figure 5. Distribution of the ages of X-ray selected PMS ob- 
jects a s inte rpolated from the isochrones derived from the Siess 
et al. (12000? ) evolutionary tracks. PMS candidates with A(R — 
Ha) 0.12 are indicated by the dashed line, whilst those with 
A(R — Ha) < 0.12 are indicated by the solid line. The total num- 
bers of objects with A(R - Ha) ^ 0.12 and A(R - Ha) < 0.12 
are respectively 70 and 302. 1-cr Poisson uncertainties have been 
indicated for both distributions. 

least 10 counts per energy channel. In the following, we dis- 
cuss only the spectra of those sources that have an aver- 
age, vignetting and exposure-corrected count rat^3 above 
5 x 10" 3 or 10~ 2 ctss _1 for EPIC-MOS or pn respectively 
and are not associated with OB stars. The X-ray spectra of 
the OB stars were discussed in a separate paper (Paper II). 
We use the numbering scheme of the sources as introduced 
in Paper I. 

For sources that display no significant variability, we 
have analyzed the merged spectra of all observations, while 
for those sources that experience a flare in their light curve, 
we rather analyze the spectra of the observation when the 
flare occurred. The spectra were then modelled by means of 
rather simple models: absorbed mekal opticall y thin ther- 
mal plasma model (Mewe et al. ll985l ; Kaastra 1 19921 ) with 
one or two temperatures, absorbed power law, absorbed 
bremsstrahlung model or absorbed blackbody energy dis- 
tribution (see Tables[2]to|l). 

Most of the flaring sources display a spectrum that is 
best represented by a single temperature mekal model with 
best fit temperatures of kT ^ 3.4 keV. For the majority 
of the non-flaring bright X-ray sources, a second temper- 
ature is needed to model the spectra. In these cases, we 
find that most of the objects (except sources #478 and 
#535) have moderate hydrogen column densities (Ah ^ 
0.36 x 10 22 cm -2 , in reasonable agreement with the aver- 
age Nh estimated from the cluster reddening, see below) 
and temperatures of order 0.7 and 2.5 - 3.5 keV for the soft 

4 We note that the actual (uncorrected) count rates are substan- 
tially lower for sources near the edge of the field of view. 
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and hard component respectively. The flaring sources thus 
appear to have significantly harder spectra than the non- 
flaring objects. Finally, some of the objects have spectra 
that are best modelled by a power law spectral energy dis- 
tribution. Three out of five sources in this group have no 
optical counterpart (not even in the deeper extension of the 
SSB06 catalogue). 

3.2 X-ray light curves 

Light curves were extracted for each source over the same 
area as for the spectra. They were background subtracted 
and corrected for the effect of the good time intervals. We 
used different energy bands and different time bins from 10 s 
to 5 ks and we applie d vari ous tests to search for variability 
(see e.g. Sana et al. |2004| ) . Considering the typical count 
rates of the sources studied here, the most useful time bin 
to investigate variability is 1 ks and we thus focus here on 
the light curves obtained for this binning. 

3.2.1 Comments on individual sources 

In this section, we discuss the light curves and spectra of 
X-ray sources that are either variable or have ambiguous 
spectral fits. Sources #181, 228, 251, 269, 290, 297, 350, 
407, 478, 535, 568 and 603 display no significant variability 
in their light curve and are not discussed here. Their spectra 
are however described in Tables[5] to U 

3.2.1.1 #41 This source displays little variability. The 
spectrum is best fitted by an absorbed power law model 
(Table©, although a blackbody model (with kT = 1.50 keV 
and TVh < 0.08 x 10 22 cm -2 ) yields a fit that is only slightly 
poorer (xl = 0.96). 

3.2.1.2 #43 The merged spectrum of this source dis- 
plays a clear Fe K line. The best quality fit is achieved for 
a single temperature mekal model (see Table©. A power 
law fit yields xl = 1-08 for iV H = 0.39 x 10 22 cm" 2 and 
r = 1.65. The source is somewhat brighter and the spec- 
trum somewhat softer during observation #6. 

3.2.1.3 #48 This object displays a flare during observa- 
tion 4. 

3.2.1.4 #100 While no X-ray emission is detected dur- 
ing observations 1, 2, 4 and 6, the source undergoes a strong 
flare during observation 3 (see Fig. [6} and is also clearly de- 
tected (though with a much lower count rate) during obser- 
vation 5. The spectra obtained during observation 3 display 
a moderate Fe K line and are best fitted with a kT ~ 4.3 keV 
mekal model (see Table©, though a power law model with 
a photon index of 2.2 also yields a xl of 0.90. The spectrum 
from observation 5 is of much lower quality, but still sug- 
gests that the source was much softer (kT ~ 0.9 keV) and 
about 50 times fainter compared to observation 3. 

3.2.1.5 #138 The X-ray flux is about 3 times larger 
during observation 5, though the light curve exhibits no clear 
flaring behaviour. 



Light curve of source #100 during obs ff3 
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Figure 6. X-ray flare of source #100 as observed during observa- 
tion 3 with EPIC-MOS 2 (top panel) and EPIC-pn (lower panel). 
The light curve was extracted for pulse invariant channels 500 to 
10000 (i.e. for energies approximately in the range 0.5 - lOkeV). 
The time is given in ks from the beginning of the first observation. 



3.2.1.6 #171 No significant variability is found for this 
source. The spectrum cannot be fitted by a single temper- 
ature mekal model, though a 2-T model (Nu = 0.79 x 
10 22 cm~ 2 , kTi = 0.27keV and kT 2 = 1.61 keV) yields 
xl = 1.30, slightly better than the power law model quoted 
in Table© 

3.2.1.7 #175 This source is clearly detected in all our 
observations. During the first three pointings, the aver- 
age MOS count rate is about 6 x 10 -4 ctss _1 . At some 
point between observation 3 and 4, the X-ray flux of the 
source has increased by at least a factor 30 (see Fig.©. 
Indeed, during the fourth pointing, the count rate slowly 
decreases from ~ 2.0 x 10~ 2 to ~ 0.6 x 10 _2 ctss _1 . Dur- 
ing the last two observations, the source is still brighter 
than at the beginning of our campaign though the aver- 
age MOS count rate decreases slowly from ~ 5.1 x 10~ 3 
to ~ 2.6 x 10 -3 ctss -1 . The spectrum of observation 4 in- 
dicates a rather hot (kT = 4.3 keV, see Table© plasma, 
though a power law fit (r ~ 2.0) is only marginally poorer 
(xl = 1.05). During observations 5 and 6, the spectrum is 
slightly softer (kT decreases to 3.4 keV or F increases to 2.6) . 



3.2.1.8 #204 A flare is seen in the MOS 1 count rate 
towards the end of observation 4. However, no pn or MOS 2 
data are available for this source to confirm this result. 

3.2.1.9 #241 The count rate of this source is about a 
factor 15 (respectively 6) larger during observation 5 (resp. 
6) than during observations 1 to 4. It is likely that the source 
has undergone a flare at some time between the end of the 
fourth and the beginning of the fifth pointing, and that our 
data actually cover the decreasing tail of this flare. 



An XMM-Newton view of NGC 6231. III. 



7 



Table 2. Best fit parameters of the absorbed single temperature mekal models (wabs*mekal) fitted to the EPIC spectra of faint X-ray 
sources in NGC 6231. The source numbering scheme is adopted from Paper I. The third and fourth columns indicate the combinations of 
observations and instruments used in the spectral analysis. Columns [9] and [10] yield the observed and absorption-corrected X-ray flux 
respectively. The fluxes are expressed in 10 — 14 erg cm - 2 s — 1 and are evaluated over the 0.5 - lO.OkeV energy band. Col. 11 yields the 
internal numbering scheme of SBL98. Numbers between brackets are rather taken from the most recent extension of the SSB06 catalogue 
(see Paper I) and mostly correspond to sources out of the SBL98 field of investigations. 



Source IDs 




Obs. 


Inst, 








kT 


xl 


d.o.f. 


f x 8 


f x" 


SBL98 


V 


V - Iq 


XMM UJ 


Paper I 






(10 


22 cm -2j 




(kcV) 
















[1] 


[2] 


[3] 


[4] 




[5] 




[6] 


[7] 


[8] 


[9] 


[10] 


[11] 


[12] 


[13] 


165334.5-415858 


#43 


all 


M2 + pn 





»t:Sl 


8 


96 + 2.22 
-1.48 


1.05 


440 


22.4 


26.6 








165335.9-413807 


#48 


#4 


M2 





»t:S? 


3 


43+2.83 
-1.20 


1.62 


10 


22.4 


26.5 


[1987] 


[15.89] 


[1.47] 


165348.3-414008 


#100 


#3 


M2 + pn 





19 +.06 
-.05 


4 


32 + 0.89 

-0.69 

26+1.08 
-0.73 


0.87 


79 


24.2 


28.6 


4664 


20.69 


2.90 


165358.6-414643 


#175 


#4 


Ml + M2 





1 ^+ 07 
-.06 


4 


1.01 


66 


37.9 


43.5 


5433 


18.27 


2.08 


165410.8-413918 


#285 


#1 


EPIC 





,o+04 

-.05 

17+-° = 
— .05 


4 


44 +0.78 
-0.55 
+ 0.62 
-0.54 


1.19 


153 


34.2 


41.1 


[3321] 


[17.36] 


[2.04] 


165429.2-415047 


#442 


#1 


EPIC 





4 


1.29 


118 


13.9 


16.2 


7554 


20.24 


2.18 


165431.9-415107 


#467 


#3 


EPIC 





,,+ .06 
-.06 

, 7 +06 
-.06 


5 


70 +1.80 
-0.95 


1.09 


129 


21.8 


26.2 


7757 


18.12 


1.95 


165432.0-415155 


#469 


#6 


EPIC 





4 


,7+0.76 
-0.59 


1.20 


105 


18.1 


22.6 


7772 


17.46 


1.87 


165456.8-420156 


#583 


all 


pn 





11± .X5 


5. 


+ 10.61 
58 -2.23 


0.89 


27 


3.2 


3.4 


[25121] 


[18.08] 


[2.21] 



Table 3. Same as Table(2] but for the absorbed power law models (wabs*power). 



Source IDs Obs. Inst. TVtt V v d.o.f. f° bfc f^ olr SBL98 V V — I n 

,, , " XX t- 

XMM UJ Paper I (10^ cm Z ) 

[1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] 



165334.4 


414424 


#41 


all 


Ml + M2 





+ .23 
-.23 





89+- 16 
° -.13 





94 


179 


19 


1 


21 


5 








165358.0 


415220 


#171 


all 


EPIC 





47+ 03 
-.03 
,.+ .07 


3 


42 +.17 

"±i!S 

-1.16 


1 


45 


284 


4 


3 


13 


7 








165404.9 


415645 


#228 


all 


Ml + pn 





3 





97 


100 


3 





7 


3 


332 


15.71 


1.49 


165405.9 


414509 


#241 


#5 


EPIC 





37+09 

-.07 

59+ 21 
-.15 


2 





76 


80 


16 


1 


24 


8 








165450.6 


420201 


#568 


all 


pn 





3 


«1« 





87 


49 


2 


7 


10 


1 


[4893] 


[16.17] 


[2.00] 



EPIC light curve of source § 175 




140 160 240 250 260 320 340 

Time elapsed [1000s] 



Figure 7. EPIC-MOS 1 (top panels) and MOS 2 (lower panels) 
light curve (for PI in the range 500 to 10000) of source #175 
during XMM observations 3, 4 and 5 (from left to right). The 
time is given in ks from the beginning of observation 1. 



3.2.1.10 #249 No significant variability is found for this 
source during the different pointings, although the EPIC 
count rates are somewhat larger during the last observation. 



3.2.1.11 #254 During observation 5, the source exhibits 
a small flare with an increase of the count rate by a factor 5 
to 6 which is consistently seen in all three EPIC instruments. 

3.2.1.12 #258 The X-ray flux is largest during obser- 
vation 1 and decreases slowly afterwards. 

3.2.1.13 #285 This source experiences a strong flare 
during the first observation and a subsequent decrease of 
the count rate during the second observation (see Fig. [8}. A 
small flare is seen in all EPIC instruments during the fifth 
pointing. The spectrum obtained during the flare of the first 
observation can be fitted by a 1-T mekal model (see Tabled, 
although a fit of equal quality can be obtained for a power 
law model with iV H = 0.37 x 10 22 cnT 2 and Y = 2.10. 

3.2.1.14 #304 A strong flare is seen during the second 
observation where the count rates increase by about a fac- 
tor 100 (Fig.[9]). The spectrum during the flare exhibits a 
prominent Fe K line. 

3.2.1.15 #324 The source is about a factor 2 brighter 
during the fourth observation, whilst the count rates of the 
different instruments during the other observations exhibit 
no coherent significant variability. 

3.2.1.16 #442 A flare with an increase of the count rate 
by a factor 10 occurs during the first observation. During 
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Table 4. Same as Table[2] but for the absorbed 2-T mekal models (wabs*(mekali+mekal2)). 



s. 
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1 
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0.87 
1.05 
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1.05 
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1.00 
0.81 
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0.86 
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62 
254 
173 

56 
214 

69 
124 
192 
334 
131 
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168 
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92 
110 
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3.6 
5.2 
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9.6 
3.0 
4.8 
6.0 
2.7 
90.1 
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3.7 
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EPIC light, curve of source # 285 
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10 20 30 60 70 80 320 340 




10 20 30 60 70 80 320 340 
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Figure 8. EPIC-MOS 1 (top panels), MOS 2 (middle panels) and 
pn (lower panels) light curve (for PI in the range 500 to 10000) of 
source #285 during observations 1, 2 and 5 (from left to right). 
The time is given in ks from the beginning of observation 1. 



the subsequent observations, the count rate slowly decreases 
(Fig.Hnj). 

3.2.1.17 #467 The source exhibits a flare during the 
third observation (Fig, lll[) . During the flare, the spectrum is 
equally well described by a kT = 5.7 keV single temperature 



EPIC light curve of source § 304 
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Figure 9. Same as Fig.[8]but for source #304 during observations 
1, 2 and 3 (from left to right). 



mekal model (see Table[2} or a power law model with V — 
1.81 and iV H = 0.34 x 10 22 cm" 2 . 

3.2.1.18 #469 Another flaring source with an increase 
of the count rate by a factor ~ 60 during the last observation 
(Fig. ll2[l . During the flare, the spectrum is again equally 
well described by a single temperature mekal model (see 
Table[2]) or a power law model with T = 2.09 and = 
0.38 x 10 22 cm -2 . 
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EPIC light curve of source § 442 



EPIC light, curve of source # 469 
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Figure 10. Same as Fig. [8] but for source #442 during observa- 
tions 1, 2 and 3 (from left to right). 



Figure 12. Same as Fig.|Ubut for source #469 during observa- 
tions 4, 5 and 6 (from left to right). 
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EPIC light curve of source § 467 
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Figure 11. Same as Fig. [8] but for source #467 during observa- 
tions 2, 3 and 4 (from left to right). 



3.2.1.19 #553 This source displays a steady increase of 
its count rate over the duration of the second observation 
(Fig. ll3[l . The light curve is truncated at the end of the ob- 
servation due to the rejection of a soft proton flare from our 
data (see Paper I). During the remaining observations, lit- 
tle variability is seen. While we adopt the 2-T mekal model 
(see Tabled}, which fits the Fe K line better, we note that a 
power law model with T = 2.27 and Nb 
yields a fit of formally equal quality. 



0.42 x 10 22 cm 



3.2.1.20 #583 A flare is seen during the fourth observa- 
tion (Fig. ll4[) . The spectrum has a rather low quality because 
of a low actual count rate. 



3.2.2 X-ray flares 



Most of the flares displayed in Figs. 161 14l exhibit an impulsive 
rise of the source X-ray luminosity within a few ks followed 
by a slower decay. The exceptions are sources #100 and 
probably #553 that display a much slower rise on timescales 
of 10 ks and ^ 20 ks respectively. Other sources (such as 
#467) show a slower fl ux inc rease than expected for typi- 
cal flares. Stelzer et al. (Il999l ) suggested that such a shape 
could be explained by a rotational modulation of the X-ray 
emission produced by the flare. In such a configuration, the 
flare actually occurs on the hidden face of the star and then 
becomes slowly visible with the rotation of the object. It is 
however difficult to investigate more this possibility in the 
present case as the rotational periods of our objects are not 
known. We finally note that the detected flaring sources dis- 
play a harder spectrum than quieter sources, with an average 
temperature during flares about 4-4.5 keV. 

Some quantitative insight into the flare events can be 
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Table 5. Properties of the flaring sources: r yields the 1/e decay time, kT b s is the best fit plasma temperature during the flare (a 
value between brackets corresponds to the hottest component of a 2-T fit) and I is the loop half-length. For comparison, the stellar radii 
computed from the effective temperatures and bolometric magnitudes derived in Sect. [2] are provided in the sixth column. Column 7 
indicates the value of the R — Ha index compared to the value of this index for main sequence stars with the same intrinsic V — Iq 
(whenever available). Finally, the last two entries correspond to the JHK colours derived from the 2MASS measurements. 



Source 


Obs. 


T 

(ks) 


(keV) 


I 

(Ro) 


R* 
(Re) 


A(R - Ha) 


J-H 


H - K 


#48 


4 


8.8 


3.43 


1.0 






0.641 


0.215 


#100 


3 


5.3 


4.32 


0.7 


1.3 








#204 


4 




(3.64) 
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0.03 


0.620 


0.171 


#254 


5 


7.3 


(3.69) 




3.7 


0.08 


0.388 


0.067 


#285 


1 


9.4 


4.44 


1.2 






0.788 


0.343 


#304 


2 


5.5 


(8.2) 




1.6 




0.204 


0.071 


#442 


1 


7.6 


4.04 


0.9 


0.8 








#467 


3 


16.0 


5.70 


2.4 


1.7 


0.45 


0.991 


0.432 


#469 


6 
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EPIC light curve of source # 553 
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Figure 13. Same as Fig. [8] but for source #553 during observa- 
tions 1, 2 and 3 (from left to right). 



obtained assuming that the flaring plasma is confined in a 
closed coronal loop as is the case for solar flares. For these 
situations, Serio et al. (|l99ll ) have established a simple ana- 
lytical relation between the loop half-length the maximum 
temperature at the top of the loop T max and the thermody- 
namic decay time r: 



ax 

where I is expressed in R©, kT in keV and r in ks. This 
relation is valid for the initial decay phase during which the 
emission measure decay is roughly exponential and has been 
applied to flares in a variety of active stars (e.g. Briggs & Pye 




240 250 260 320 
Time elapsed [1000s] 



340 



Figure 14. EPIC-pn light curve (for PI in the range 500 to 10000) 
of source #583 during observations 3, 4 and 5 (from left to right). 
The time is given in ks from the beginning of observation 1. 



120031 . Giardino et al. l2004l and Favata l2005l ). When the tem- 
peratures are determined from spectral fits of EPIC data, a 
correction has to be applied to the observed temperature to 
derive the max imum temperature a t the top of the loop (see 
Giardino et al. |2004| , Briggs & Pve l2003r ). 

Whenever possible, we have determined the 1/e decay 
times of the flares by a least square fit of an exponential 
decay to the decay phase of the observed light curves. Usu- 
ally, the results for the different EPIC instruments are in 
good agreement (to better than 10% of r), except for sources 
#254, #285, #442 for which we have a larger dispersion of 
~ 20% between the values determined from the three light 
curves. For sources #175 and #241, our data do not cover 
the maximum of the flare and are hence not suited to con- 
strain the decay time, although we note that the light curve 
of #241 during observations 5 and 6 is well fitted with a sin- 
gle exponential decay suggesting a rather large decay time 
30 ks). 

Table[5]yields the properties of the flaring X-ray sources 
with a sufficiently large count rate to study their spectrum 
and light curve. All flaring sources have rather large tem- 
peratures during the flare. For those spectra that were fitted 
with a 2-T model, Table[S] indicates the temperature of the 
hotter model component between brackets. The decay times 
range from about 2 ks to 16 ks and the loop half-lengths are 
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in the range 0.3 to 2.4 Rq. Most loops appear rather com- 
pact (Z ^ R t ) exce pt perhaps for source #467. 

Favata l]2005ft notes that while most flaring structures 
are indeed relatively compact, extremely slowly decaying 
flares have been detected in some young stellar objects 
(with r ~ 40 ks) whi ch co rrespond to rather large loops 
(I ~ 14 Rq). Favata (|2005l 'l accordingly suggests that the 
magnetic structures that produce these large loops link the 
star to the circumstellar disk. We note that, although our 
dataset is not suited to clearly identify slowly decaying 
flares, the decay seen for source #241 might be a candi- 
date for such an event. Unfortunately, the source could not 
be associated with an optical counterpart though a near- 
IR counterpart was found. Hence, we can not exclude that 
#241 is a background source. 

In this context, it is also interesting to note that among 
the five objects in Table[5]for which a A(_R — Ha) value can 
be computed, only source #467 shows a strong Ha emis- 
sion. Further evidence for circumstellar material around this 
source comes from its near-IR excess. From Table[5j it ap- 
pears that #467 is also the source with the largest loop 
half-length (I ~ 1.4 R.*) in our sample. 



3.3 X-ray fluxes 

Using the bolometric magnitudes derived above from the 
SBL98 photometry, and the vignetting and background cor- 
rected EPIC count rates determined by the source detec- 
tion algorithm, we can build an empirical Lx/iboi relation 
for the sources in the NGC 6231 field. In order to convert 
the EPIC count rates into X-ray luminosities in the 0.5 - 
lOkeV range, we took the average of the available count 
rates of the pn, MOS1 and MOS2 instruments assuming a 
conversion factor 2 between the MOS and pn count rates. 
Next, we computed a flux/pn count rate conversion factor 
assuming that the sources have an absorbed optically thin 
thermal plasma spectrum with a temperature of 1 keV and 
JV h = 5 8 x 10 21 x E(B - V) = 2.7 x 10 21 cm" 2 (Bohlin et 
al. Il978h . Note that adopting a temperature of 2keV (more 
typical of the spectra of the brightest sources) would only 
slightly increase this conversion factor (by 14%) . Finally, we 
converted the unabsorbed fluxes into luminosities, assuming 
that all sources are members of NGC 6231 (DM = 11.0). 
Figure [T5l displays the resulting logix/I/boi as a function of 
Mboi. This value increases from about —7 for BO stars (the 
hottest objects plotted in this figure) to almost —2 for the 
lower mass objects. A si milar trend was already observed in 
NGC 6530 (Rauw et al. 12002ft . 

Of course the results in Fig. [15] could be somewhat bi- 
ased by the fact that some of the sources are detected only 
during their flaring activity. We have repeated the same pro- 
cedure for the brighter sources discussed hereabove. In the 
bottom panel of Fig. 1151 we distinguish between the flaring 
and non flaring sources. We find that the average luminosity 
during a flare can reach up to 3 x 10 32 ergs -1 (source #304) 
and logLx/iboi can temporarily increase to values as large 
as —1.00 (source #100). On the other hand, the brighter 
steady sources also display rather large X-ray luminosities, 
up to ~ 4 x 10 31 ergs -1 , as well as large log Lx/iboi, up to 
~ —2.4. Therefore, we conclude that even non-flaring PMS 
sources in NGC 6231 can be rather X-ray bright. 
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Figure 15. Top panel: logarithm of Lx/^bol as a function of 
the bolometric magnitude for the EPIC sources with optical coun- 
terparts in the SBL98 catalogue (see text). Lx refers to the X-ray 
luminosity in the energy range 0.5 to 10 keV. The various symbols 
correspond to different combinations of EPIC instruments used 
in the determination of the X-ray luminosity (see the insert in the 
upper left corner). Bottom panel: log Lx/^bol as a function of 
the bolometric magnitude for those EPIC sources with optical 
counterparts in the SBL98 catalogue and for which a detailed 
modelling of the spectra was possible (see Tables [21 EH and [4j . 
Asterisks correspond to sources that have a rather constant flux 
during our six observations, whereas hexagons stand for flaring 



4 DISCUSSION 

Our XMM-Newton campaign provides the deepest X-ray 
view ever of NGC 6231. While the X-ray images of the clus- 
ter are dominated by its high-mass members (see Paper I), 
we detect a substantial number of sources that have faint 
optical counterparts. The typical spectra of these sources, 
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their flaring behaviour, their X-ray luminosities as well as 
the location of their optical counterparts in the H-R diagram 
indicate that these sources form a population of PMS stars. 
Compared to pre vious optical investigations of the cluster 
(Sung et al. 1 19981 ). the number of PMS stars is increased by 
a factor 20 at least. 

An interesting feature is the age spread among the low- 
mass PMS objects (~ llMyr) that exceeds the main se- 
quence lifetime of the massive O stars in the cluster. This 
featur e was already noted by previous authors. Raboud et al. 
l|l997h obtained photometry of NGC 6231 in the Geneva sys- 
tem. Fitting isochrones to the main sequence and post-main 
sequence stars believed to be single, these authors derived 
an age between 3.2 and 4.5 Myr for the more massive stars, 
whilst they noted that the PMS stars detected in their obser- 
vations were mostly older than 3.2 Myr and had ages up to 
10 Myr. A similar situation was found by Sung et al. (1998) 
who derived an age of 2.5 - 4 Myr for the massive stars in 
the cluster, while they found a n age spread of 11 Myr in the 
PMS objects. Raboud et al. (|l997l ) accordingly suggested 
that there was a continuous low-mass star formation activ- 
ity in NGC 6231 at previous times, which came to an end 
about 3.2 - 4.5 Myr ago when the massive O stars of the clus- 
ter formed during a nearly coeval burst. From our results, 
the age distribution of X-ray selected low-mass stars peaks 
around 1 - 4 Myr, i.e. at the time of the formation of the 
early-type cluster members. However, the formation of low 
mass stars in NGC 6231 started already more than 10 Myr 
ago at a rather slow rate. The star formation rate subse- 
quently increased, culminating in a starburst event during 
which most of the PMS as well as the more massive objects 
formed. Then, about 1 Myr ago, the formation of stars in 
NGC 6231 came to an end rather abruptly. 

In the case of NGC 6530, Damiani et al. (|2004T ) report 
on a spatial age gradient, suggesting a sequence of star for- 
mation events. To check whether a similar situation holds 
for NGC 6231, we have plotted the spatial distribution of 
X-ray selected PMS stars of different ages in Fig. 1161 It can 
be seen that, in the case of NGC 6231, there exists no clear 
separation between the different groups. While this could be 
a projection effect (if star formation were progressing along 
our line of sight towards the cluster), we nevertheless con- 
clude that there is no evidence in our data for a spatial age 
gradient in NGC 6231. 

The star formation history of NGC 6231 provides con- 
straints on future modelling of massive star formation in 
open clusters. Several properties of the NGC 6231 cluster 
could be explained by a massive star formation scenario in- 
volving physical collisions in a dense cluster core (Bonnell 
et al. I199ST ). These features are the large binary fraction 
among the massive stars in this cluster, the spatial concen- 
tration of PMS stars around the massive cluster members 
and last but not least, the age difference between the PMS 
objects and the more massive stars. Indeed, N-body simu- 
lations of compact cluster cores including the effects of stel- 
lar ev olution and physical collisions (Portegies Zwart et al. 
1 19991 ) show that, because of mass segregation, physical colli- 
sions between stars are frequent in the cluster core. The re- 
sulting mergers could lead to the formation of very massive 
(^ 100 Mq) stars within 3-4 Myr. The stars are rejuvenated 
by each new collision, therefore considerably extending their 
lifetime. The products of such a process would thus appear 
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Figure 16. Spatial distribution of the X-ray selected low- 
mass stars in NGC 6231. The triangles, squares, stars and filled 
hexagons indicate stars with ages of less than 2 Myr, between 2 
and 5 Myr, between 5 and 9 Myr and more than 9 Myr respec- 
tively. 

as some sort of 'blue stragglers' in the H-R diagram of the 
cluster. 

Alternatively, a moderate age difference could also re- 
sult if massive star formation proc eeds t hrough accretion 
as propo sed b y Norberg & Maeder (|2000l ) and Behrend & 
Maeder i|200ll) . These authors assumed that the PMS ob- 
jects accrete a constant fraction of the observed outflow 
rates, which in turn increases with the stellar luminosity 
and hence the stellar mass. In their models, the stars be- 
come visible on the birthline once a significant part of the 
surrounding cocoon has been disper sed. For instance, in the 
models of Norberg & Maeder ( 2000), an accreting protostar 
that will eventually form a 120 Mq star would spent half of 
the accretion time below a mass of about 2Mq. The time 
needed for an accreting protostar to increase its mass from 2 
to 120 Mq would then be about 1 Myr. These models require 
rather high accretion rates (of order 8 x_10 -3 Mq yr _1 to 
form an 85 M star, Behrend & Maeder l200lh . Though this 
value seems high, it is appar ently not forbidden by the re- 
sults of Wol fire fc Cassinelli (119871 ). The models of Behrend 
& Maeder (|200ll ) accordingly predict that low mass stars 
form first and high mass stars form slightly later. However, 
the duration for the formation of an 80 Mq star is only about 
10% larger than that of an 8 Mq star. Henceforth, this sce- 
nario alone cannot explain the large age spread observed in 
Fig.0 



5 CONCLUSIONS 

In this third paper of the series, we discussed the properties 
of the NGC 6231 X-ray sources with faint optical counter- 
parts. Their typical X-ray spectra and luminosities, their 
flaring behaviour as well as their location in the H-R dia- 
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gram indicate that the bulk of these objects are most prob- 
ably PMS stars with masses between 0.3 and 3.0 Mq. Their 
observed age distribution suggests that the low-mass stars 
in NGC 6231 started their formation at a relatively slow rate 
over 10 Myr ago. This process then culminated 1 to 4 Myr 
ago in a starburst-like event that coincided with the time of 
formation of the massive stars in the cluster. No difference 
in the location of the cTTs and wTTs in the cluster as well 
as no spatial age gradient were observed. 

The X-ray spectra of the brightest PMS sources are 
well reproduced by rather hard thermal plasma models. 
About one third of them displays flaring activity during the 
XMM-Newton campaign. The observed logLx/Lboi values 
increase with decreasing bolometric luminosity up to values 
of logLx/iboi ~ —2.4 for non-flaring sources and to even 
larger values during flares. 

Finally, we qualitatively discussed possible implications 
for the massive star formation in NGC 6231. Briefly a sce- 
nario implying the merging of low-mass objects to create 
high-mass stars would probably allow to explain most of 
the observational properties of the cluster. In particular, it 
will likely be able to reproduce the large binary fraction seen 
in the massive star population as well as the important time 
(several Myr) separating the beginning of formation of the 
low- and high-mass stars. 

To conclude, we outline that X-ray emission is probably 
one of the best selection criteria for PMS objects. Beside a 
detailed understanding of the massive star population, we 
now have probably identified a significant fraction of the 
PMS population of the cluster. The determination of the 
initial mass function of both the high-mass and low-mass 
star populations could bring further constraints on the star 
formation in NGC 6231 and, more generally, on the possible 
link between the formation of these two populations. This 
task is deferred to a forthcoming paper in the series. 
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